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We present the first large sample of accurate iron abundances and temperatures for 

clusters at redshifts > 0.14. We find that the Fe abundance shows little or no evolution 
Li^ ■ out to z ~ 0.3. This and the early formation epoch of elliptical galaxies in clusters 

[ — ■ indicate that most of the enrichment of the intracluster medium occurred at z > 1. 

If clusters represent fair samples of the universe, then global metal production is 2-5 
times greater than is inferred from recent studies of galaxies at low and high redshifts. 



Subject headings: galaxies: abundances, galaxies: formation, intergalactic medium 



1. Introduction 

The hot intergalactic medium in low redshift clusters of galaxies has been enriched in iron 
to ~ 1/3 solar (Edge &: Stewart 1991, Yamashita 1992) by processes that, if understood, are a 
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powerful diagnostic of galaxy formation. Since the mass in the x-ray emitting gas exceeds that in 
the stars in galaxies by factors of 2-10 (assuming H a = 50 km s _1 Mpc -1 as we do throughout 
this paper; we also assume £1 = 1), the IGM is the main reservoir of metals in clusters. Moreover, 
the strong correlation of iron mass with optical light ( Arnaud et al. 1992 ) implies a coupling of 
the constraints on the evolution of metals from the x-ray data on high redshift clusters and the 
optical data on the evolution of elliptical galaxies. If, as many theorists believe, clusters are fair 
samples of the universe, as well as the only systems likely to have retained all the metals produced 
by stars, then studies of the evolution of the IGM metallicity strongly constrain the star formation 
history in the universe. 

There were few measurements of Fe in clusters at z > 0.1 before the launch of the ASCA 
satellite ([Tanaka et al. 1994|). Neither early results from ASCA (|Mushotzky 1996, Fsuru et 



al. 1996| ), nor subsequent detailed individual measurements ( Bautz et al. 1994 , Donahue 1996 



Matsuura et al. 1996J ), showed evolution in the Fe abundance out to z ~ 0.3. Here we report on a 



large sample of high signal-to-noise ASCA spectra of clusters at z > 0.14, and find no evidence for 
evolution of the cluster metal abundance. Since recent optical data reveals only passive evolution 
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in cluster ellipticals out to z ~ 0.4, most of the cluster metals were produced at z > 1. We thus 
infer much higher star formation rates in the z > 1 universe than derived from optical/UV surveys. 



2. Observations 

We extracted and analyzed ~ 40 spectra of clusters at z > 0.14 from data in the ASCA 
public archives before 12/96. We concentrate on the 21 high signal-to-noise observations with 
90% confidence error in Fe abundance < 0.12 - a moderately large sample upon which simple 
and unambiguous statistical tests can be made. At kT ~ 6 keV, the product of the flux, f x , 
and exposure time must exceed ~ 1.4 10 -7 erg cm -2 to attain this accuracy. For a typical 
ASCA exposure time of 40 ks (net exposures for two clusters in our sample, MS1358 and AC118, 
were ~ 70 ks), a minimum flux f m in = 3.5 10 -12 erg cm -2 s _1 is thus required. Before the 
release of the ROSAT all-sky survey, such objects were rare: e.g., in the IPC database only 19 
clusters have z > 0.14 and f x > f m in- Cross-correlation of the RASS and Abell catalogs yields 
50 clusters with f x > 1.4/ m j n at z > 0.14 ( [Ebeling et al. 1996| ), but further cross-correlation 
with the ASCA observation log identifies only 19 clusters with f x > f m in- Thus we include most 
of the high signal-to-noise, high-z clusters. (The flux from one cluster that apparently met our 
criterion, Abell 1774, turned out to be dominated by an AGN.) This subsample (Table 1) spans 
the bolometric x-ray luminosity (about twice the 2-10 keV luminosity) range from 10 45,1 to 10 46 ' 3 
erg s , overlapping the low-z sample measured by previous satellites. 



We applied standard processing to the data (e.g., Mushotzky et al. 1996] ); and, for both 



the SIS and GIS, used a 6' extraction radius, corresponding to 3 Mpc at z ~ 0.3 - sufficiently 
extended to include the total flux. At z ~ 0.15, 10-30% of the flux may fall outside the extraction 
radius. We used all-sky backgrounds accumulated during the ASCA deep survey pointings for 
clusters with z < 0.2, and background from the same field for the higher z systems. Spectra were 
fit to an isothermal, collisional equilibrium plasma model at the optical redshift of the source. The 
redshifts obtained from ASCA data are within 0.013 of the optical values (e.g., for the relatively 
faint MS1358 the fitted redshift is z fit = 0.323(+. 011,-0.022) compared to z opt = 0.327, while for 
the brighter ZW3146 z fit = 0.292 (+0.004,-0.0061) and z opt = 0.291). 

Cooling-flow-plus-isothermal-plasma models were used in cases (ZW3146, A2390) where a 
cooling flow is clearly present in the ROSAT image (cf. |AUen et al. 1996] ). As found by [Allen et~aL 



1996, the extra component has little effect on the derived abundances of the isothermal plasma. 
Our fits are in excellent agreement with results for those clusters in Table 1 that were previously 
analyzed: RXJ 1347 flSchindler et al. 1996Q , ZW3146, A1835 and E1455 QAllen et al. T996J ) , A1204 



( Matsuura et al. 1996[) , and MS0451 ( ponahue 1996J ), (The results are new for all others.) The 



parameters derived from ASCA spectra (Figure 1) are robust for these high flux objects. 



Table 1. High Signal-to-Noise Clusters 



Cluster log L x & log kT (keV) Z Ff ° 8Z_ 

A1204 45.08 0.587 

A1722 45.23 0.768 

A1246 45.18 0.798 

MS2137-24 45.25 0.643 

A963 45.30 0.829 

MS1358+62 45.30 0.813 

A1413 45.36 0.797 

A1763 45.38 0.953 

A2218 45.41 0.857 

A773 45.45 0.952 

MS1455 45.47 0.660 

A520 45.48 0.934 

MS0451-03 45.60 1.01 

AC118 45.65 0.970 

A2390 45.73 0.949 

ZW3146 45.74 0.802 

A1689 45.77 0.947 

A2204 45.78 0.903 

A2219 45.92 1.07 

A1835 45.92 0.910 

RXJ1347-11 46.28 0.968 

a Log bolometric x-ray luminosity in erg s^ 1 assuming q = and H — 50 km s _1 Mpc -1 . 
b Fe abundance with respect to the Anders & Grevesse 1989 photospheric abundances. 
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c The error in the Fe abundance - 90% confidence for one parameter (Ax = 2.7) 
d Optically determined redshift. 



3. Cluster Metallicity as a Function of Redshift 

Our sample of 21 clusters has nearly identical Fe abundance properties to the low-z Ginga 
sample ( Yamashita 1992] , Butcher 1995), i.e. a mean abundance (Fe) = 0.27 ± 0.013 and 



variance 0.076 (Figure 2) compared to (Fe) = 0.29 ± 0.014 and variance 0.062 (all errors are 90% 
confidence for one parameter of interest). A K-S test shows no significant difference between the 
samples. We find no evidence of systematic abundance variations with luminosity. The weak 
abundance-temperature correlation (Figure 3) is not significant at 90% confidence. We conclude 
that this sample, with (z) = 0.26, shows no statistical differences in Fe abundance from a similar 
sample of low-z objects. Results for the systems with larger abundance errors agree with the 
above conclusions. 

Although our sample is not flux- or redshift-limited, it has no obvious selection biases and 
should be representative of z ~ 0.25 clusters. Its mean luminosity, (log Lf, i) = 45.4, and mean 
temperature, (kT) = 7.5 keV, are slightly higher than in the low-z sample drawn from a flux 
limited survey ((log Lt, i) = 45.16, (kT) = 5.96 keV). The temperature-luminosity relationships of 
the z ~ 0.25 and low-z samples show no significant differences flMushotzky fe Scharf 1997 ). 



Discussion 



In high-S7 hierarchical clustering scenarios ( Kauffmann 1995 ), the mass of rich clusters grows 



by ~ 2 from z ~ 0.2 to the present. Moreover, systems at z ~ 0.4 evolve differently from a cluster 
of the same richness observed today, with high-z clusters assembled over a shorter time interval 
and undergoing more merging in the few Gyr prior to the epoch when they are observed. Thus 
the present and low-z samples should differ if the accreted material has distinct properties, and 
the surprising constancy of both the mean and variance of the Fe distribution over this redshift 
range (corresponding to a lookback time of ~l/4-l/3 the age of the universe) argues for a very 
homogeneous accretion process or a low-S7 universe where cluster growth is essentially completed 
at z ~ fl^ 1 ( White 1996| ). The substantial spread in cluster metallicity at z = and z = 0.25 is 



unexpected within purely gravity-driven hierarchical clustering scenarios. Many of these clusters 
include an excess of blue galaxies compared to low-z systems (the "Butcher-Oemler effect"). None 
of these have anomalous metallicity: whatever processes converted high-z blue spiral galaxies into 
low-z spheroidal systems had no major impact on the IGM. 

The absence of evolution over the past ~ 5 Gyr implies early enrichment of the IGM. Thus, 
ram pressure stripping of infalling galaxies at low-z cannot dominate the enrichment, but much of 
the cluster Fe could be produced by SNIa ( Renzini et al. 19"93| ): reconciling large Fe- mass-to-light 



ratios with modest present-day SNIa rates requires the rate to have a short rise time followed 
by steep negative evolution such that 90% of the Fe is produced at z > 0.5. Most of the cluster 
metals, perhaps including most of the Fe, were produced at high z by SNII ( [Loewenstein & 
Mushotzky T996D. 



The stars in elliptical galaxies that we see today come from the same population that produced 



the metals in the IGM (e.g., Elbaz et al. 1995); moreover, the elliptical galaxies in several rich 



z ~ 0.3-0.4 clusters of galaxies show only passive evolution (e.g., Bender et al. 199q ). Therefore, 
the star formation in elliptical galaxies that produced the cluster metals must have stopped > 2 
Gyr before z ~ 0.4, placing the period of cluster metal formation at z > 1 — in excellent agreement 
with the lack of evolution inferred from ASCA data. This is in contrast to the history of metal 
production by disk galaxies, where star formation peaked at 1 < z < 2 and continues through the 
present epoch flPei fc Fall 1995| , pVladau et al. 19961 ). 



The containment by the deep cluster potential of the metals synthesized and ejected by 
supernovae during the epoch when the stars that make up present-day elliptical galaxies formed, 
enables us to calculate the z > 1 metal production by (proto-)ellipticals and estimate their 
contribution to the total metals produced at z > 1 and over the Hubble time. Assuming that 
relative abundances in clusters are constant, the mass in metals produced per unit (present-day) 
blue-band elliptical galaxy luminosity is ~ 0.3 in solar units. If the star formation (though not 
necessarily the dynamical evolution) in elliptical galaxy precursors proceeded in a similar way 
in the field and in clusters, then the product of the metal mass-to-light ratio and the early-type 



galaxy luminosity density ( Loveday et al. 1992 ) yields the total mass density of metals from 
early- type galaxies, 

p z = 1.4 W 7 M Q Mpc- 3 . (1) 



This is consistent with the value, derived by|Zepf fc Silk 1996, required to account for the increase 



in stellar mass-to-light ratio with luminosity as due to stellar remnants from primordial star 
formation with an IMF biased towards high mass stars. 

The metals ejected by proto-ellipticals, as calculated above, exceed what is estimated to be 



locked inside the stars in normal massive galaxies today by 2-5 ( Madau et al. 1996 ). Since the 
latter is consistent with the star formation rate derived from galactic UV light, star formation in 
proto-ellipticals must have proceeded in an obscuring environment or in sub-luminous fragments 



( Loewenstein &: Mushotzky 1996| ). Since half the star formation in disks has occurred since 



2 



1 but, star formation in ellipticals was completed by z = 1 (cf. |Cowie et al. 1996J ), 



proto-spheroids increasingly dominate star formation at early epochs. If elliptical galaxy star 
formation predominantly occurred at 1 < z < 6, the average metal production rate was 

^ = 3.4 10" 3 Moyr^Mpc- 3 , (2) 

at 

compared to averages from high-z, UV surveys of galaxies of 0.89, 0.46, and 0.17 xlO~ 3 M Q yr _1 
Mpc~ 3 for z > 1, z > 2, and z > 3, respectively. 

Thus, global star formation rates derived from UV surveys are severe underestimates. At 



z > 0.3, where Ha redshifts out of the ground-based optical band, [OII]3727 or the UV continuum 
are used as star formation indicators. Recent ASTRO-2 UIT and HUT observations ( [Fanelli et 



al. 1996| ) show that there is often a shift of a factor 2 in specific star formation rates using the 



UV continuum instead of Ha; and furthermore, comparison of IRAS 60/i luminosities with Ha 
values indicates another factor of 2 offset (see also, Buat 199^ , Deharveng et al. 1994). Since even 



small amounts of reddening are sufficient to remove objects from the high-z galaxy samples and to 
change their inferred star formation rates, we believe that the cluster data are much more robust 
indicators of the true metal formation history of the universe. 



5. Conclusions 

The Fe abundance of the intergalactic medium in rich clusters of galaxies undergoes little 
if any evolution from z ~ 0.3 to z ~ 0. This strongly constrains high-fi hierarchical theories in 
which clusters rapidly grow during this period. A real and constant spread in the metallicity at 
all observed redshifts implies stochasticity, not only in the construction of these massive systems, 
but also in the efficiency and/or timing of star formation within their constituent galaxies. 

A self-consistent history of star formation in disk galaxies has emerged from a synthesis of 
observations of damped La systems, deep ground based and Hubble Space Telescope surveys of star 
forming galaxies, and spectral synthesis and chemical evolution modeling. However, this census of 
star formation in the universe is incomplete, since most of the metals produced at high z were 
expelled from galaxies. A complementary picture of star formation in proto-spheroids is emerging, 
driven by the lack of evolution in IGM metallicity reported here and optical evidence that together 
indicate that cluster metals originate from the precursors of elliptical galaxies at z > 1. The 
L-dependence of M/L ( Zepf fe Silk 199(f ) and gross and relative IGM elemental abundances 



( Loewenstein fe Mushotzky 199q ) are explained if this star formation produces an IMF skewed 



toward high masses, precipitating a massive outflow of metal-enriched gas. This star formation 
must occur in sub-luminous fragments or be dust-enshrouded to avoid detection by UV surveys, 
but the IR and/or UV background produced is within reach of near-future measurements (Zcpi 



& Silk 1996). The IGM of rich clusters are reservoirs for primordial metal-rich galactic outflows 
driven by this star formation, and their investigation with ASCA has enabled us to discover that 
~ 4 times more metals were created at z > 1 than are inferred from high redshift galaxy data. 

We thank Prof. Y. Tanaka for early encouragement of this work, and acknowledge the efforts 
of the entire ASCA team. 
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Fig. 1. — 68, 90 and 99% error contours for the temperature and Fe abundance in three rich clusters, 
as indicated, illustrating clusters with small, typical, and large abundance uncertainties. 

Fig. 2. — Fe abundance vs redshift. The error bars show 90% confidence uncertainty limits. 

Fig. 3. — Fe abundance vs gas temperature 
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